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ROCKET-POWERED MODEL INVESTIGATION OF LIFT, DRAG, AND 
STABILITY OF A BODY-TAIL CONFIGURATION AT 
‘MACH NUMBERS FROM 0.8 ТО 2.5 AND 
ANGLES OF ATTACK BETWEEN #6.59 


By Warren Gillespie, Jr., and Albert Е. Dietz 
SUMMARY 


A configuration having a body of fineness ratio 16.6, an unswept 
horizontal teil of aspect ratio 2.7, and a highly swept vertical tail 
was &eropulsed continuously in pitch during free flight with and with- 
out а sustainer rocket motor operating. The Mach number range covered 
by the test was 0.8 to 2.5 and the model angle of attack did not exceed 
и. £6.59. Zero-lift-drag and drag-due-to-lift data were obtained during 
л coasting flight of the model. Normal-force, pitching-moment, and static- 

stability data were obtained with and without the rocket motor thrusting. 


Correlation of power-off normal-force- and pitching-moment-curve 
slopes with theoretical estimates based on slender-body theory as pre- 
sented іп МАСА RM A52D29 gives good agreement with experiment. The 
reciprocal of the normal-force-curve slope closely approximates experi- 
mental drag-due-to-lift values. The amplitude of the angle-of-attack 
oscillations during flight with power on was approximately 50 percent 
smaller than the amplitude during coasting flight, and the normal-force- 
and pitching-moment-curve slopes were effectively higher during powered 
flight. The above differences between coasting and powered-flight 
values may be largely due to cross-coupling effects induced by opera- 
tion of the sustainer rocket motor and the asymmetry of the cruciform 
tail. 


INTRODUCTION 


Ав part of a general program to determine the aerodynamic charac- 
қ teristics of wing-body-tail combinations at supersonic speeds and 
moderate angles of attack, a body-tail configuration was flight-tested 
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at the Langley Pilotless Aircraft Research Station, Wallops Island, Va. 
The aeropulse technique as presented in reference 1 and a large solid- 
fuel sustainer rocket motor were used to obtain lift, drag, and static 
stability of the model at Mach numbers from 0.8 to 2.3 and angles of 
attack up to +6.5°. The model consisted of a body of fineness ratio 16.6 
with an unswept horizontal pulsing tail of.aspect ratio 2.7 and a highly 
swept fixed vertical tail of aspect ratio 1.08. The horizontal tail was 
aerodynamically actuated between stop settings of +2.0° in approximately 
а square-wave pattern with a frequency of from 2 to 5 cycles a second. 
The basic aerodynamic parameters in pitch were determined from the 
response of the model to the tail motion. 


SYMBOLS 
CN normal-force coefficient, Dont 
Ce chord-force coefficient, - 2L 
Ср drag coefficient, Са cos a + Cy sina 
Cm pitehing-moment coefficient, Nt - 

27.5954 

an normal acceleration, ft/sec? 
a7 longitudinal acceleration, ft/sec? 
g acceleration due to gravity, 32.2 ft/sec? 
а dynamic pressure, lb/sq ft 
M Mach number 
R Reynolds number, where the reference length is 1 ft 
S body cross-sectional area, 0.267 sq ft 
d body diameter, 0.585 ft 
о, angle of attack, deg 
Ө angle of pitch, deg А | 0l с 
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Y angle of yaw, deg 

p roll rate, radians/sec 

б all-movable horizontal-tail deflection, #2.09 
€ effective upwash at the tail, deg 

Iy moment of inertia about Y-axis, slug-ft© 

Iy moment of inertia about X-axis, slug-ft? 

Iz moment of inertia about Z-axis, slug-ft* 

P rocket chamber pressure, lb/sq in. abs 

T rocket motor thrust, lb 

t time, sec 

x distance from center of pressure to model base, in. 
1 model length, 116.15 in. 

m model mass, slug 

Subseripts: 

g ground 

f flight 


Note: One апа two đots over symbols депође first and second time 
derivatives, respectively. 


MODEL 


А two-view drawing of the test configuration is shown in figure 1. 
The fuselage was a body of revolution of fineness ratio 16.6. Ordinates 
defining the nose shape are given in table I. The geometric and mass 
cheracteristics are listed in table II, and photographs are shown in 
figure 2. The maximum body-diameter-tail-span ratio was 0.250 for the 
vertical tail and 0.315 for the horizontal tail. The horizontal tail 
was mass-balanced and pivoted £2.09 about the 0.55-exposed-mean- 
aerodynamic-chord point. 


n eu X ONETDENTIAT, NACA RM L5hCOh 


The model was of metal construction and carried & solid-fuel 
sustainer rocket motor in addition to ап eight-channel telemeter with 
angle-of-attack, pressure, апа accelerometer instruments. 


The model and its booster are pictured in the launching attitude 
in figure 5. Total impulse was approximately 17,000 10-вес for the 
solid-fuel Deacon booster and 7,660 lb-sec for the sustainer motor. ` 


TEST 


Data were obtained after the model separated from the booster. 
During flight of the model alone, a square-wave pulse was continuously 
generated as the tail automatically flipped between stop settings due 
to a reversal in direction of the tail lift. 


The quantities measured by the telemeter system were normal and 
longitudinal accelerations, angle of attack, horizontal-tail position, 
free-stream total pressure, and sustainer-rocket-chamber pressure. The 
angle-of-attack indicator mounted on the vertical tail became inopera- 
tive before sufficient data could be recorded. Тһе velocity obtained 
from Doppler radar was used in conjunction with tracking radar and 
radiosonde data to calculate Mach number, Reynolds number, and dynamic 
pressure experienced by the model during flight. The variation of 
Reynolds number per foot length and dynamic pressure with Mach number 
are shown in figure 4. The model experienced а coasting period before 
and after the period of flight with the sustalner motor thrusting. 


Aerodynamic coefficients were determined for coasting and power-on 
flight conditions. Data at а particular Mach number were obtained only 
during constant tail deflections of £2.09. Coefficients are based on 
the body maximum cross-sectional area of 0.267 square foot and on the 
body diameter of 0.58 foot. A detailed explanation of data analysis 
is given in reference 1. Pitching moment was determined from the dif- 
ference in normal accelerations measured at the nose and near the center 
of gravity of the model. 


ACCURACY 


The random error in the data is indicated by the scatter of the 
experimental points which is generally much less at the highest Mach 
numbers. The maximum absolute accuracy of a quantity obtained from а, 
single instrument is usually better than 2 percent of the total cali- 
brated instrument range. The probable error is approximately 1 percent. 
Presented on the following расе are the ranges of the telemeter instru- 
ments used in the test model: БЕТІНЕН 
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Nose angle-of-sttack indicator, deg . ............... 18 
Normal accelerometer at the nose, g-units . ............ £8 
Normal accelerometer near the center of gravity, eine V Me ж SLO 
Longitudinal accelerometer, coasting, g-units . . . . . . . . . 1 to -T 
Longitudinal accelerometer, power on, g-units’. . . . . . . . -5 to 60 
Rocket chamber pressure, Ib/sq in. . . so . ......... О to 1400 
Free-stream total pressure, Ib/sq in.  . .. ........ -5 to 115 


An additional source of inaccuracy in the final results may be cross- 
coupling effects of induced rolling апа yswing motion. This possibility 
is discussed more thoroughly in & later section. 


RESULTS AND DISCUSSION 


Model pitching response.- Typical portions of the pitching response 
of the model are shown in figures 5, 6, and T. Figure 5 shows а part of 
the decelerating Plight time history of the first coasting period and 
indicates how the model angle of attack, normal acceleration, and Mach 
number varied as the tail moved in an approximate square-wave pattern 
of £2.0° deflection. Figure 6 shows portions of the angle-of-attack 
telemeter record throughout the flight of the model after separation 
‘from the booster. The aeropulse motion developed immediately after the 
model became free of the booster. The oscillations during the first 
coast period were uniformly symmetrical. Тһе record shows that the 
angle-of-attack indicator at the tail started to vibrate at Mach num- 
ber 1.2 and then ceased to function. The record also shows that the 
amplitude of the sinusoidal oscillations was suddenly reduced and varied 
after the sustainer motor started thrusting. Figure 7 shows more clearly 
that the amplitude of the angle-of-attack oscillations was reduced 
approximately 50 percent during operation of the sustainer motor. After 
sustainer burnout the oscillations increased in amplitude but were 
asymmetrical during the early part of the second coasting period and 
irregular during the latter part. 


Cross coupling.- The curious variations that occurred in the angle- 
of-attack response of the test model may be largely due to cross-coupling 
effects of rolling and yawing motions. Such motions may have been 
instigated by operation of the sustainer rocket motor and amplified by 
the asymmetry of the horizontal-tail and vertical-tail plan forms and 
horizontal-tail deflections of £2.00. At the time of sustainer ignition, 
it is estimated that thrust misalinement in yaw could have started a 
yawing motion with an initial maximum yaw angle of approximately two 
times the amount of angular thrust offset. А thrust offset of approxi- 
mately O.h? in pitch is estimated from the trim offset indicated by 
figure 7. From reference 2 it may be determined that the value of the 
rolling parameter C1g/9 for the horizontal tail of the test model is 
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small compared with that for the vertical tail. The test model there- 
fore may have experienced an aerodynamic rolling moment at combined 
angles of pitch and yaw. This is different from the case of a symmetri- 
cal cruciform missile for which the rolling tendencies of the individual 
tail surfaces would cancel at angles of attack up to about 6°, as shown 
by the experimental data of reference 3. 


An investigation of the effect of steady rolling on longitudinal 
and directional stability was reported in reference 4 where the stability 
was studied by means of the Euler moment equations. Іп particular the 
equation for pitching moment M = түб - (Iv, - Ty )ўр shows the term 
introduced by the combination of rolling and yawing velocity. This 
moment term may have affected the pitching moment measured by the 
accelerometers of the test model. The accelerometer measuring normal 
force may have been similarly affected by a term containing the product 
of yaw angle and roll velocity as shown by the following expression 


7. = mV(& ~ 6) - шур 


Calculations were therefore made to determine the possible effect 
of cross coupling on the transient angle-of-attack response of the test 
model. The method of calculation is given in appendix A. The results 
of these calculations are shown in figure 8 and indicate that for reason- 
able values of roll rate and initial yawing velocity cross coupling 
probably caused the asymmetrical pitching oscillations of the second 
coasting period. Cross coupling may also have been largely responsible 
for reducing the amplitude of oscillation during powered flight. Addi- 
tional calculations made without cross-coupling terms show that during 
the second coasting period when the tail setting was ~2.0° the experi- 
mental oscillation can be very exactly matched by the calculation method 
of reference 1. During part of the second coasting period, therefore, 
the effect of cross coupling on the pitching motion was apparently 
negligible. 


Sustainer motor.- Operation of the sustainer rocket motor in flight 
was not appreciably affected by the high longitudinal and normal flight 
accelerations encountered. This is shown in figure 9 where the rocket 
chamber pressure obtained during a static ground. test is compared with 
the chamber pressure measured in flight. Thrust during flight was 
calculated from the expression Tp = Т„(Ре/Р;) where the values used 


occurred at corresponding burning times during the ground and flight 
tests. Thrust corrections to account for differences in back pressure 
during the ground and flight tests were negligible. 


The effect of the thrust lift component T sina on the amplitude 


of the pitching response during powered flight was also considered. 
The method of calculation is given in appendix B. Тһе calculation 
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shows that the thrust lift component probably reduced the power-on 
pitching amplitude approximately ЈЕ percent below the amplitude without 
thrust acting. 


Angle of attack at zero normal force and at tail flip.- Figure 10 


shows the angle of attack at the instant the model pitched through zero 
normal force during coasting and powered flight. The dashed line is 
interpolated for a tail setting of 09 and lies close to zerc angle of 
attack in the low supersonic range. As the Mach number approached 2.0, 
the angle of attack for Cy = О changed approximately 0.59. The angle- 
of-attack indicator had a Zon delta surface. It is probable that the 
indicator changed trim as the flow velocity normal to the leading edge 
of the indicator surface changed from subsonic to supersonic velocity 
at Mach number 2.0. The trim effect of small pressure disturbances 
arising from minute imperfections in the construction of the indicator 
would probably be different, depending on whether the leading edge was 
subsonic or supersonic. 


The angle of attack at the instant that the horizontal tail started 
to flip from one stop setting to the other is presented in figure 11. 
Assuming that at this instant the effective angle of attack at the tail 
was zero, the effective upwash at the tail can be approximately estimated 
from the relation : 


€ = -Gflip - 5 


The average effective upwash over the tail obtained from the two tail 
settings was 1.19 at an average angle of attack of 0.99 over the Mach 
number range of the tests. This is greater than would be expected from 
two-dimensional potential-flow theory and the experimental data of fig- 
ure 20 of reference 5. The high upwash indicated by the tail flipping 
may be due to the degree of mass balance of the tail which could affect 
the start of the tail flipping. 


Drag.- A typical drag polar is presented in figure 12. The polar 
is composed of data from one oscillation with tail settings of +2.0°. 
The drag characteristics are summarized in figure 15 which shows the 
variation of zero-lift drag and drag due to Lift with Mach number. The 
zero-lift drag coefficient reached a peak value òf 0.52 at Mach num- 
ber 1.0 for this configuration and uniformly decreased at higher Mach 
numbers. The drag due to lift had a minimum value of 0.046 at Mach 
number 1.0 and increased with increasing Mach number. The param- 


eter TE is in good. agreement with drag-due-to-1ift values. 
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Normal force, pitching moment, and tail effectiveness.- Figures 15 
to 18 present typical plots of normal-force and pitching-moment coefficients 
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end summarize the variation of the normal-force- and pitching-moment- 
curve slopes with Mach number. The unusual variation of the normal 

force and pitching moment with angle of attack exhibited in figure 15 

for a teil setting of 2.0° occurred during the asymmetrical oscillations 
of the second coasting period. The change іп slope of these curves 
between points B to D may be due to cross-coupling effects of induced 
rolling and yawing motions during this portion of the oscillation. 

Slopes obtained between points D to E and F to A were used in prepara- 
tion of figures 16 вий l7. Figures 16 and 17 show that the normal- 
Ғогсе- and pitching-moment-curve slopes due to angle of attack were 
effectively higher during powered flight than during coasting flight. 

The increase in these slope values between powered. and coasting flight 

is of the same magnitude as that indicated by regions B to D of figure 15. 
It appears therefore that the cross-coupling terms mVyp and (Ig - Ту){р 
discussed previously may have caused the differences between coasting and. 
powered-flight values. 


Figure 16 also shows the tail lift effectiveness. No difference 
between power-on and power-off values of this parameter can be determined 
because of the rather large scatter of the data. 


Figure 18 presents a comparison of normal-force-curve slope and 
tail effectiveness with the theory of reference 6. Very good agreement 
is noted for the normal-force-curve slope at all Mach numbers and the 
tall effectiveness at the highest Mach numbers of the test. 


Model aerodynamic center.- Figure 19 shows the variation of model 
aerodynamic center with Mach number for both power-off and power-on 
flight conditions. There is very little difference due to power effects. 
The most rearward power-off location of the aerodynamic center at 21 рег- 
cent body length from the base occurred near Mach number 1.0. At higher 
Mach numbers the aerodynamic center moved uniformly forward as the tail 
effectiveness decreased. Figure 20 compares aerodynamic-center location 
calculated by the method of reference 6 with the experimental curve of 
the present test. Using the first-order theory of reference 7 to deter- 
mine the nose contribution rather than slender-body theory gives very 
good agreement between the curves. 


Induced jet -effect.- The possibility of induced jet effects on the 
powered-flight values of Су, Ome and aerodynamic-center location 


was also considered. However, reference 8 shows that induced jet effects 
on the afterbody of a rocket model would not be expected to occur above 

а Mach number of 1.4. Below Mach number 1.4 strong jet effects would 
probably not extend more than 1 base diameter forward of the jet exit. 
The horizontal tail of the test model was located over 1 base diameter 
ahead of the jet exit. Furthermore, if induced jet effects at the model 
tail had caused the Cy, and Ст, slope differences shown in figures 16 
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and 17, then the model aerodynamic center should have shifted farther 
rearward than the small amount shown in figure 19. It seems probable 
that large induced jet effects were not present. 


CONCLUSIONS 


An investigation of lift, drag, and stability of a rocket-propelled 
body-tail configuration between Mach numbers of 0.8 and 2.5 with and 
without the sustainer motor operating leads to the following conclusions: 


1. During coasting Plight, the drag due to lift increased with 
increasing Mach number throughout the supersonic range of the test and 
was slightly higher than the reciprocal of the normal~force-curve slope. 
The lift and aerodynamic-center location were accurately predicted by 
use of the interference factors derived from slender-body theory. 


2. The amplitude of the pitching oscillations was smaller during 
powered flight and the normal-force- and pitching-moment-curve slopes 
were effectively higher than values obtained during coasting flight. 
These differences may be largely due to cross-coupling effects induced 
by operation of the sustainer rocket motor and the asymmetry of the 
eruciform tail. 


Langley Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., February 18, 1951. 
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APPENDIX A 
EFFECT OF CROSS COUPLING ON TRANSIENT ANGLE OF ATTACK 


The calculation of the effect of cross coupling on the transient 
angle-of-attack response to a step input tail deflection of the test 
model was made in the following manner. In order to simplify the calcu- 
lations, only angular motions in pitch and yaw were considered.  Aero- 
dynamic damping was neglected. The aerodynamic moments for equal angles 
of pitch and yaw were assumed equal. The model inertia about the longi- 
tudinal axis was assumed zero and roll rate was assumed constant. The 
moment equations in pitch and yaw with axes fixed in the model reduced 
to 


E NEA А 
57.5 a 37.5 py Cm Фава Cy. 598 


15, кы 17. • 
TUI КЕ a DNE NAE — И ай Sd = 
57.5 Y * 57,5 қайысы 


or, Since Тү = ір for the test model, 


Е 
z 
Я 
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= 
+ 
d 
= 


and nim 7 ~ 5 and р = = Тһе solution with no damping is 


2C - 20za,° + 20 
a(t) = аут + jm i сов 21. 5а1% + 
Бај“ - 50681 2 + Ст 


20-842 ~ 2a оС 
"nce dM y sin 5T. 2845 + ng - + 205 сов 57. aot + 
584 - 30624? + Сч Ба - Den + С? 


202827 - 20182 


X 


oe ааа DT. Jagt 
бас” ~ 5Сбар + CT 
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where 


6 C 
ap? = 26 + <8) - С" 


2 
Сі = а(0) 
Co = а(0) 


“г (о) - 2 2 Е 
б = рҰ(0) - да(0) + ра(0) + = 


о, = Ру у(о) - 2 а(о) 


= С D 
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APPENDIX B 


EFFECT OF THRUST LIFT COMPONENT ON MAXIMJM 


AEROPULSE ANGLE OF ATTACK 


The effect of the thrust lift component T sina on the maximum 


aeropulse angle-of-attack response was calculated using the logarithmic 
chart of figure 21. Тһе chart was originally derived by Mr. Robert і. ` 


Nelson of the NACA Langley Pilotless Aircraft Research Division and is 
based on the equations of longitudinal motion with freedom in pitch 
about the y-axis and translation in the z-direction. In the derivation 
of the chart the assumptions of instantaneous tail flipping and equal 
but opposite slopes at succeeding times, 50; t4, . . ., of tall flip 
are made. | 


The solution for the maximum aeropulse response is 


ab | 
a = 1 ~ | - шау, WE віп 57.2аЉу + cos б.а) - 


_ атыр) 25) ~b(tmtti) in БТ. 
L TE l+ SR e gin 57.5аб sin 57.58%1 
-bt : 
1 + е *(cos DI, Jaty - 5. sin 57.3at1) 


and 


) 
аза - Zb + Сн ста) 


where 


тү а 
_ 57.398 ` Ош ay 25 
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са E 
Ош By Na + 27.546 "a 
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-bt 
e lsin „Зад, 
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1 + E сан 51.25а54 


Note that tj is the time at the tail flip following the previous maxi- 
mum amplitude at the time ty. 
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TABLE І.- CONTOUR ORDINATES OF NOSE 


Station, in. Body radius, 
from nose in. 


95 
1 
2 
6 
7 
Т 
О 
22 
2 
2 
2 


ХАН Ун го о оого 
~ сом ном 


оо со 


16 ~ , 
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TABLE ІІ.- CHARACTERISTICS OF MODEL 


Воду: 
Maximum diameter, Ғ9........ 
Base diameter, ft ......... 
Length, ft . . . . . © © © e we we on n 
Body fineness ratio ........ 
Nose fineness ratio . .... 0.» 
Boattail fineness ru SG UE ae E 
Boattail angle, deg . . . . ex « . 


Horizontal tail: 
Span, total, ft .......... 
Exposed area, sq ft ........ 
Aspect ratio, total span and area . 
Aspect ratio, exposed span and area 
Taper ratio, exposed . . .. « « « « 
Sweepback at 0.5с, deg . . a . e sœ 
Airfoil section . . . « e « в © e e 


Vertical tail: 
opan, total, ҒСГ.......... 
Chord at center line . . . « « « e « 
Aspect ratio, total span and area . 
Sweepback at L.E., deg . . ea . . . 
Sweepback at T.E., deg . . . . e e e 
Airfoil section . $9 e s w e «a e a è 


Model weight, Jb: 
With sustainer rocket loaded... 
With sustainer rocket empty . .. 


* 


Moment of inerti& in pitch or уам, Blug-ft?: 


With sustainer rocket losded . . .. 
With sustalner rocket empty » e s. à 


1/4-inch beveled 


Center of gravity with sustainer rocket loaded or empty, 


percent length from base . . . . e à 


омм Бо о CH 
У оол л о ЕЛ 


сло © бого © 


: hepercent hexagonal 


flat plate 


183.6 
144.9 


Sta Sta Sta 
6600 8913 102.88 


§ 

о 

99 Model c.g. 

Ps power on and off 
| 


Sta 
Gr) @ 


· Коонбт ин VOVN 


Ё 

$ 

® 
500 diam 


Totat-pressure tube 056 = : 
© Side view Sta Sta 15° 
8990 99,38 


Figure 1.- Two-view drawing of the test configuration. ALL dimensions 
are in inches. 
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(a) Side view. 


1.-85555 


(ъ) Bottom view. 


Figure 2.- Photographs of test configuration. 
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Figure 5.- Model and booster on launcher. 
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(а) Reynolds number. 


Reynolds number , R 
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Figure 3 .- Variation of Reynolds number and dynamic pressure with 
Mach number. 
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Figure 5.- A typical portion of the flight time history. 
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(a) Development of sinusoidal oscillation following model separation from booster. 
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(b) Flutter of indicator at tail; М-1.24, 


Figure 6.- Typical angle-of-attack oscillations during flight of model. 
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(c) Reduction in amplitude following ignition of sustainer motor at 72 seconds 
flight time. 
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(d) Asymmetrical oscillations of increased amplitude following burnout at 8.4 
seconds flight time 
Figure 6.- Continued. 
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(e) Transition to irregular oscillations of varying amplitude. 


Figure 6.- Concluded. 
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Figure T.- Effect of powered flight on model angle of attack. 
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Figure 8.~ Effect of constant roll and initial yaw rates on angle-of- 
attack response to tail deflection of 2.09, М = 2.02. 
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Figure 9.- Comparison of thrust developed during ground and flight testes. 
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Figure 11.- Model angle of attack at tail flipping. 


хооцст ИЧ VOVN 


бё 


е МАСА RM 154С0Щ. 


5.2 


-8 ` О ôs - 200 deg 
D 8: 2.00 deg 


МАС 
А RM ІБІКСОВ 


51. 


wi 
th Mach number 


Figur 
e 15 
„- Vari 
lation 
of zer 
o-lift drag 
coeffi 
cient 
ала 
dra, 
g due 
to li 
ft 


22 


Wer 


TERMED 
УДА | T 
COPE srr, 


NACA RM L5hCOh 


Е 
| 
Р 
a ЋЕ 
Жү 
5 An 
; 
p wd 
Е gf 
з 59 
ЈЕ 
1 
Бъ 
4% 
Н 
d 


МАСА ЕМ L5hCOL тнщ 


el Lab 1° 
zl 9, 


HP EHH, 
ГТ Pee 
70, 
0" 
ТҮГЕ, 
к) — yu m i 


См 0 


(b) Power-on flight. 


Figure 14.- Continued. 
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Figure 14.- Concluded. 
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(с) Second coasting period. M = 1.90. 
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(а) Normal-force coefficient against angle of attack. 
Figure 15.- Typical variation of normal-force and pitching-moment 


coefficients during asymmetrical oscillations of second coasting 
period. 
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(b) Pitching-moment coefficient against angle of attack. 


Figure 15.- Concluded. 
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Figure 17.- Variation of pitching-moment-curve slope with Mach number. 
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Figure 18.- Comparison of theory and experiment for Cy, and Ср 
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Figure 20.- Comparison of theory and experiment for aerodynamic-center 
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